Nonequilibrium current driven by a step voltage pulse: an exact solution 
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One of the most important problems in nanoelectronic device theory is to estimate how fast or how 
slow a quantum device can turn on/off a current. For an arbitrary noninteracting phase-coherent 
device scattering region connected to the outside world by leads, we have derived an exact solution 
for the nonequilibrium, nonlinear, and time-dependent current driven by both up- and down-step 
pulsed voltages. Our analysis is based on the Keldysh nonequilibrium Green's functions formalism 
where the electronic structure of the leads as well as the scattering region are treated on an equal 
footing. A model calculation for a quantum dot with a Lorentzian linewidth function shows that 
the time-dependent current dynamics display interesting finite-bandwidth effects not captured by 
the commonly used wideband approximation. 
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Understanding coherent transport of charge and spin 
through the scattering region of a nanoscale device is 
the central problem of nanoelectronics theory which has 
tremendous scientific and technological importance 0. 
From both theory and application points of view, an im- 
portant issue which has yet to be resolved is to predict 
how fast or how slow a nanoelectronic device can turn 
on/off a current from quantum mechanical first princi- 
ples. Indeed, one cannot develop an electronic technol- 
ogy unless the operational speed of the device can be 
designed and controlled. This issue is closely related to 
the transient transport phenomenon, which is becoming 
an extremely important problem of nanoelectronic device 
physics, as can be observed in such effects as photon- 
assisted tunneling 3, electron turnstiles ^ and ringing 
behavior in the time-dependent current [j, la. Recent 
real-time measurements of electron dynamics kf\ have fur- 
ther raised interest for the study of transient quantum 
transport. 

The purpose of this work is to investigate transient 
quantum transport far from equilibrium, for nanoelec- 
tronic systems in the Lead-Device-Lead (LDL) config- 
uration where "Device" indicates the scattering region 
which is connected to the outside world by the leads. 
When the LDL system is driven far from equilibrium by 
a step-shaped voltage pulse, we discovered an exact so- 
lution to the time-dependent current J{t), thereby an 
exact transient quantum transport picture is obtained. 
To the best of our knowledge, this is the first time an 
exact solution is found for far-from-equilibrium transient 
transport dynamics in the quantum regime, and it pro- 
vides a valuable and unambiguous physical picture of how 
charge current is turned on and off by a bias voltage pulse 
through devices of the LDL form. Importantly, the exact 
solution provides, for example, the correct current decay 
time scale after the bias voltage is turned off, and the 
transient current that follows the turning-on of a con- 
stant voltage. 



For quantum devices in the form of LDL, the the- 
oretical formalism best suited to the study of time- 
dependent transport is, perhaps, the Keldysh nonequilib- 
rium Green's functions In the NEGF formalism, 
the time-dependent current through the phase-coherent 
scattering region of the device is given in terms of lo- 
cal Green's functions. However, these Green's functions 
cannot, in general, be solved analytically. Previous stud- 
ies of time-dependent transport have thus so far relied 
on the so-called wideband limit (WBL) which is a 
simplifying assumption where the coupling between the 
scattering region and the external leads is taken to be in- 
dependent of energy. In other words, the WBL neglects 
any electronic structure of the device leads. Therefore, 
when the electronic structure of the leads is important, 
i.e. for leads with finite bandwidth such as those made 
of semiconductors, nanotubes, nanowires, etc., a theory 
beyond the WBL is necessary. In this regard, numerical 
approaches have been put forward in Refs. |^ 

Our starting point is the formalism described in Refs. 
0, 0, The Hamiltonian of the LDL device is 

a mn 

(tka,nci^dn + tka,n4cka) , (1) 
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where cj^.^^ (cka) with a — L,R creates (destroys) an 
electron with momentum k in the left (L) or right (i?) 
lead, and (dn) creates (destroys) an electron in a 
single-particle state labeled by n in the scattering re- 
gion. Quantities ika.n describe coupling of the leads to 
the scattering region of the device. Chemical potentials 
in both leads are set to zero. As in Refs. 0, S fiol 
when an external time-dependent voltage is applied to 
drive a current through the device, we assume that the 
single-particle energies acquire a time-dependent shift: 



eka(t) 



Aa{t) and em„(t) = e^„ + A„„(i). It has 
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been shown [5| that the charge current through lead a is 
given by 

Ut) = -2ej' dt' I ^ImTr{e-(*-*')e*/>'^^°(*^) 

xr,(e)[G<(i,i') + /(e)G^(i,0]}, (2) 

where /(e) = {e^'^ + 1)^^ is the Fermi function, 
the hnewidth function Ta{e) has matrix elements 
(e) = 27rpa(e)t* „(e)tQ^„(e) where pa{e) is the 
density of states in lead a, the lesser Green's function 
G<{t,t') has matrix elements G<^{t,t') = i{d\^{t')dy,{t)) 
and the retarded Green's function G^{t,t') has ma- 
trix elements G^^{t,t') = -ie{t - t'){{dr,{t),dl{t')}). 
The retarded Green's function is given by the so- 
lution of the Dyson equation [H [13 G^{t,t') = 
G(^{t,t') + Jdtijdt2 G^(i,ti)S«(ti,t2)G^(i2,t'), where 
G^{t,t') is the Green's function of the scattering 
region without any leads, while the Keldysh equa- 
tion yields the lesser Green's function G'^{t,t') — 
J dh J dt2G^{t,ti)^<{ti,t2)G'^{t2,t'), where G^(i,t') 
is the advanced Green's function and S^'<(ti,i2) are 
the retarded and lesser self-energies, respectively. From 
a mathematical point of view, the main effect of bro- 
ken time-translational invariance due to the presence of 
time-dependent external fields is that the double inte- 
gral Dyson equation is not a simple Fourier convolution 
product as is the case in equilibrium or for steady-state 
transport, and therefore cannot be reduced by a Fourier 
transformation to a simple algebraic matrix equation. In 
the WBL, after one neglects the energy dependence of 
the linewidth function, the self-energy t2) becomes 

proportional to a delta function 6{ti—t2), hence G^{t, t') 
can be solved afterward 0, 0, ■ 

To solve the problem exactly without relying on the 
WBL, we define 



dt' 



From the Keldysh equation and the expression for the 
lesser self-energy [3, 

one can straightforwardly show that Eq. ^ takes the 
form 

/d 
^ImTr{r„(e)[«'„(e,i) + /(e)A„(e,i)]}, 

(4) 

where we have defined 



Therefore, once Aa{e,t) can be found, the current is en- 
tirely determined. 

To solve for Aa{e, t), we make use of an alternate ver- 
sion of the Dyson equation, namely that obtained by 
choosing only the time-dependent part of the Hamilto- 
nian as a perturbation, the unperturbed part being 
the time-independent terms. This way, the unperturbed 
Hamiltonian describes the LDL device at equilibrium, 
hence its physics is given by a time-translationally invari- 
ant Green's function G^{t — t'). A standard expansion 
of the contour ^-matrix along the Schwinger-Keldysh 
contour C H^ . followed by analytic continuation 
of the contour-ordered Green's function Gnm{T,T') = 
—i{Tc{dn{T)dl^{T')}) to the real time axis according to 
the Langreth rules yields the following alternate 

form of the Dyson equation. 



G^it,t')^G^it-t') + J dtiG^it ~ h)A{h)G^{h,t') 

+ fdh f dt2G"{t - h)V''{h,t2)G''{t2,t'), (5) 



where we define a retarded potential P V^{ti,t2) = 
[exp (-i di' A^(t')) - l] - t2) and Ef is 

the equilibrium retarded self-energy due to lead (3. Us- 
ing Eq. in the Dyson equation we obtain the 
following integral equation for the quantity Aa{e,t): 

A^ie,t) = A^ie,t) + J dt' e'<'~''^e' ^"^''^G^'it - t') 
xA{t')A^{e, t') + JdtiJ dt2 e*^(*-*^) (6) 
^e'^>'''''^''^G^(t~t,)V''it,,t2)AUe,t2). 

Here ia(e,i) = /^^ dt' e'^(*-*')e' '^'^ ^°^*^^G«(t - i') 
describes the equilibrium state and is therefore known. 
To solve Eq. we need to specify the external fields 
A{t) and (t). Here we investigate a step function pulse 
applied at time t = 0, of the form A(Q-)(t) = A(Q)6'(±t) 
where A(„) is a constant amplitude and the plus (minus) 
sign corresponds to an upward (downward) step. 

Downward step pulse. This is the situation where a 
constant bias voltage with value A^ is sharply turned off 
at time t = and remains off for subsequent times. For 
this case and from its definition, V^{ti, t2) vanishes when 
ti and t2 are simultaneously greater than zero, as well as 
when ti < t2 from the retarded self-energy. Equation © 
then takes the form 
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We see from the limits of integration that while Aa{e,t) 
is required for t > on the left-hand side of Eq. Q, 
only Aa{e, i < 0) is involved in the integrals on the right- 
hand side. This is an example of a Wiener-Hopf equa- 
tion For our LDL device under this downward step 
pulse, Aa{e,t < 0) is actually known: for t < the sys- 
tem is in steady-state under a constant bias Aq,. Hence 
Aa{e,t < 0) is easily obtained from the known steady- 
state NEGF G^(t-t') = G^{t < 0,t' < 0). Equation © 
is therefore not an integral equation but an explicit ex- 
pression for Aa{e, t) in terms of known quantities. Once 
the integrations are carried out by tedious but elemen- 
tary algebra, one obtains the following exact expression: 



C/„(e,cr) EE /:t_^{e^('+^°)*y^(t)} which are aU known 
quantities. Finally, a detailed but straightforward calcu- 
lation gives 



dz 



A,) 



2m 



A„ 



Aa 

z 



A-^A;3tf^(6,e + z)|G«(e) 



(9) 



A^{e,t) = G"\e)+ / — 



duj 



27ri uj 



^0^ 



A„ 



^0^ 



A-^A^tf^(c^,6) G«(6 + A„) 



(8) 



where we define Tf^(w,w') = [Sf (w) - -I- A^ - 

A;3)]/(w — ui' — Aa + Afs). Equation ((SJ is the first im- 
portant result of this paper. With the explicit solution 
for Aa{e,t), the time-dependent current Jait) can be ob- 
tained without further difficulty from Eq. I^J. 
Upward step pulse. In this case, Eq. © takes the 
form 

^„(e, t) = A^(e, t) + f dt' e*(^+A=)(t~t')Gi?(i _ ^'^ 



xAAc,{eA')+ f dti / ' dt2e*(^+^")(*-*i) 
Jo Jo 

xG^(t-ti)e'(^+^°)(*i-*^)T/^(ii -t2)Aa(e,i2), 

where A'^{e,t) is a known function that involves only 
Aa{e,t) and Aa{e,t < 0). This is a Volterra equation 
which has the form of a Laplace convolution product. It 
can thus be converted into an algebraic matrix equation 
by a Laplace transformation, so that Aa{e,t) is given by 
the following Bromwich integral: 



Aaie,t) 



2ni 



7 + 



dae° 



J— too 



X (1 - F„(e, a)[A + f/„(e, a)])-' a). 



where a is the Laplace variable, and we define A'^{e, a) = 
Ct^^{A'^ie,t)}, F^ie,a) = £t^,{e'(^+^°)*G«(i)} and 



where we define the contour Cj : 00-1-27 —oo+ij in the 
complex z-plane. Equation Q is the second important 
result of this work. 

Examples. The above exact solutions, Eqs. JHEJ, are 
valid for arbitrary noninteracting phase-coherent devices 
in the LDL configuration. As an concrete example, we 
now apply these results to a quantum dot (QD) with a 
single energy level eo connected to two leads described 
by a Lorentzian linewidth function Ta{e) = T'^W'^ /{e^ + 
W'^), where r° is a constant linewidth amplitude and W 
is the bandwidth. After the integrals in Eqs. H8I9|) are 
carried out by residue integration, the time-dependent 
current Ja{t) of Eq. (@J is then given by a single integral 
over all frequencies. This last integral can be easily done 
numerically as it contains /(e) which has known poles at 
the fermionic Matsubara frequencies iujn = i{2n+ 1)^/(3. 

In Fig. ^ we show the time-dependent current Jhit) 
through the left lead in response to a downward step 
pulse. At t = 0, energies in the left lead are lowered 
by Al = lOr and the energy of the QD level is low- 
ered by A = 5r where F is the total linewidth ampli- 
tude. For W » r, Al,A, the WBL result is essentially 
correct. When W becomes comparable to other energy 
scales of the problem, the WBL is seen to be a poor ap- 
proximation. The initial current Jl(0) decreases as the 
linewidth function gets narrower since less states in the 
leads are available for transport. More importantly, the 
time-dependent current can increase following the bias 
turnoff (curves (iv) and (v) in Fig. an interesting 
nonclassical behavior not displayed by the WBL current. 
A current increase after the bias is turned off was also 
observed in a previous numerical study 0. Most im- 
portantly, for devices with smaller bandwidth, the WBL 
and the exact solution predict very different time scales 
of the current decay. The inset shows an interesting os- 
cillatory behavior, not captured by the WBL, for a nar- 
row band W = 0.25r misaligned with the resonant level 
eo = -0.3r. 
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In Fig. 12 we show the time-dependent current Jl (t) 
through the left lead in response to an upward step pulse. 
In this case, the single-particle energies eta and Eq are 
suddenly raised at t = 0. Here again, the WBL is seen 
to be accurate for W ^ F, A^, A. However, interesting 
finite-bandwidth effects appear as bands in the leads get 
narrower. First of all, the asymptotic i — > cxd current 
decreases with W, corresponding to the initial current 
decrease in the downward step situation discussed in the 
last paragraph. In addition, it is seen that a positive 
voltage pulse can drive an instantaneously negative cur- 
rent, as has been observed in the numerical work of Ref. 
0. For W < A, Ai, the current oscillates around a zero 
value: the pulse drives the resonant level eo outside the 
band so that little current can flow through. 

In summary, we have presented an exact solution for 
the time-dependent current through an arbitrary non- 
interacting phase-coherent device scattering region con- 
nected to external leads with arbitrary energy-dependent 
linewidth functions, in the physically relevant case of an 
upward or downward step function voltage pulse. The 
results are general and are valid for far from equilib- 
rium transport situations. For a single-level QD with 
Lorentzian linewidth function, the WBL was seen to be 
a crude approximation to the exact solution in the case 
of narrow bands. Due to the finite-bandwidth effects, 
a number of nonclassical transient behaviors were found 
including a current increase following a downward pulse, 
a negative current driven by a positive upward pulse, and 
a vanishing asymptotic current for a finite positive bias 
pulse. The significance of our solution is the exactness of 
the results which give unambiguous nonequilibrium tran- 
sient quantum transport dynamics. From a practical ap- 
plication point of view, our formalism gives the transient 
current in terms of steady-state NEGF that can be cal- 
culated by any technique used for steady-state transport 
such as those atomistic first principles techniques of Ref. 

u 
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FIG. 2: Time-dependent current Jhit) through left lead in 
response to an upward step pulse for different band widths: 
dashed line, WBL {W ^ oo); (i) W = 20r, (ii) W = lOF, 
(iii) W = 5r, (iv) W = 2.5r, and (v) W ^ T. Units and 
parameters are the same as in Fig. Q Here again, the W = 
lOOr curve (not shown) is indistinguishible from the = oo 
curve. 



